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Superradiance from BEC vortices: a numerical study 
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The scattering of sound wave perturbations from vortex excitations of Bose-Einstein condensates 
(BEC) is investigated by numerical integration of the associated Klein-Gordon equation. It is 
found that, at sufficiently high angular speeds, sound wave-packets can extract a sizeable fraction 
of the vortex energy through a mechanism of superradiant scattering. It is conjectured that this 
superradiant regime may be detectable in BEC experiments. 
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Recent years have witnessed a growing interest in pur¬ 
suing analogue models of gravitational physics in con¬ 
densed matter systems. The rationale for such models 
traces back to a seminal observation by Unruh IJ, who 
noted a close analogy between sound wave propagation 
in an inhomogeneous background flow and field propa¬ 
gation in curved space-time. The analogy goes on by 
observing that, much like superfluid hydrodynamics is 
a large-scale effective theory of microscopic superfluids, 
field theory on a curved space-time might also be re¬ 
garded as a large-scale limit of a possible microscopic for¬ 
mulation of quantum gravity. The crucial point is that, 
whereas microscopic theories of quantum gravity are still 
largely a matter of speculation, the microscopic theory of 
superfluids is well developed. It can thus be hoped that 
the wide body of knowledge available for the latter can 
be brought to the benefit of the former fjj. For instance, 
assessing the mechanisms of sound radiation from ‘terres¬ 
trial black holes’ beyond the hydrodynamic picture may 
in principle offer new insights into the microscopic origin 
of cosmic black hole radiance, the Hawking effect, and 
other cosmological phenomena. 

A key step along this long-term program is the study 
of scattering and radiance phenomena from black holes 
whose background space-time can be associated with 
fluid excitations such as vortices. A model of fluid flow 
which seems particularly well suited to pursue the ’ana¬ 
logue gravity’ program is the so-called draining-bathtub 
geometry Q, namely a three-dimensional flow with a sink 
(vortex) at the origin. The flow field induced by the vor¬ 
tex is associated with an acoustic metric with two crucial 
ingredients of black hole physics: an event horizon and an 
ergosphere. The former is a spatial surface which allows 
only one-way propagation of physical signals (from the 
outside into the vortex), while the latter is a region from 
which part of the the vortex energy can be extracted via 
the mechanism of superradiance. 

Such a phenomenon was first studied by Zel’dovich 4) 
with regard to the generation of waves by a rotating body 
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and was then analysed as stimulated emission in black- 
hole radiance @, |]fl. Superresonance is an acoustic- 
wave version of the Penrose |8j] process, whereby a plane- 
wave solution of a scalar massless field in the black hole 
background is scattered from the ergosphere with an am¬ 
plification at the expenses of the rotational energy of the 
black hole. Such process has been shown to occur in 
a certain class of analogue (2+1 (-dimensional rotating 
black holes 0]. Later studies B, [Hi have discussed the 
frequency dependence of the amplification factor in su¬ 
perresonant scattering of acoustic perturbations from a 
rotating acoustic black hole by deriving the reflection co¬ 
efficient as a function of the frequency <Jo of the incom¬ 
ing monochromatic wave. It is found that in the range 
0 < (Jo < tnfl the reflection coefficient is greater than 
unity, with m being the azimuthal wavenumber and H 
the angular frequency of the acoustic horizon. 

The main purpose of this paper is to present a quanti¬ 
tative investigation of superradiant scattering from sonic 
holes associated with BEC-like vortex configurations. Al¬ 
though superradiant scattering from hydrodynamic vor¬ 
tices has been discussed in the recent literature mm, 
we believe that this is the first quantitative assessment of 
such a phenomenon under specific BEC-like conditions. 

In the limit of zero temperature, gaseous Bose-Einstein 
condensates are well described by the Gross-Pitaevskii 
equation (GPE) 
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where $(r, t) is the wavefunction of the condensate nor¬ 
malised to the total number of bosons N, a s being 
the s-wave scattering length and M the mass of the 
atoms. If we now use the Madelung representation 
<F(r, t) = p(r, t^ e lMe (, r d)/ h Jj]j in Eq. |QJ, where 
p(r,t) = |$(r,t)| 2 is the condensate density, the GPE 
equation takes a hydrodynamic form: the imaginary part 
is a continuity equation for an irrotational fluid flow of 
velocity v(r,t ) = V0(r, t) and density p(r, f), and the 
real part is a Hamilton-Jacobi equation whose gradient 
leads to the Euler equation. As is well known, the GPE 
is equivalent to irrotational inviscid hydrodynamics El- 
Low-frequency perturbations around the stationary 
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state are essentially sound waves (zero sound) and 
obey the Bogoliubov set of differential equations for 
the density perturbation p ^ and the phase perturba¬ 
tion in terms of the local speed of sound c(r) = 
^4:Trh 2 a s p(r)/M 2 . These equations, within the limit of 
validity of the hydrodynamic approximation can be re¬ 
duced to a single second-order equation for the phase 
perturbation [ 3 . This differential equation for 0^ = d' 
has the form of a relativistic Klein-Gordon equation 
d fl (v'-ggi^dA’) = 0, with g = det g^ in a curved 
space-time whose metric g is determined by the local 
speed of sound c and the background stationary velocity 
v. 

It should be noted that the linearization suppresses 
the quantum nature of the GPE so that, within the lin¬ 
ear perturbation theory, the circulation of vortices is not 
quantised as in BEC systems. The calculations reported 
below are aimed at examining what fraction of the en¬ 
ergy can be extracted through superradiant scattering of 
a sound wave from a vortex described by BEC-like pa¬ 
rameters. As discussed further below, the full nonlinear 
GPE will have to be used for a quantitative assessment 
of the extraction of energy quanta from BEC vortices. 

For a single vortex with a drain at r = 0 and angular 
velocity 12 in the draining-bathtub model, the velocity 
field of the flow is 

v = V0(r, </>) = car + fla 2 ^j /r. (2) 

where f and </> denote unit vectors in polar coordinates, 
a is the radius of the event horizon, and the background 
density po of the fluid and the speed of sound c are taken 
as constant throughout the flow. The acoustic metric 
associated with this configuration is 

ds 2 = — (c 2 — ((a 2 c 2 + a 4 I2 2 )/r 2 )) dt 2 + (2ca/r) dtdr 

—2f la 2 dt d(f> + dr 2 + r 2 d(j) 2 + dz 2 . (3) 

It is readily checked that this metric has an ergosphere 
whose radius is r erg = a^/\ + f l 2 a 2 /c 2 . The growth of 
the ergosphere with increasing 12 allows an increasing ex¬ 
traction of energy from the vortex in superradiance con¬ 
ditions. 

Linear perturbations of the velocity potential T satisfy 
the massless Klein-Gordon scalar wave equation on this 
background, i.e. 
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Equation 0 is solved by means of numerical meth¬ 
ods developed in [l5| for the integration of massless 
scalar-field perturbations on a rotating Kerr black-hole 


background. To this purpose, Eq. 0 is conveniently 
recast into a system of first-order (strongly) hyper¬ 
bolic equations through the definition of two conjugate 
fields S i = d'b/dx 1 and II = — {1/a) (d'S’/dt— 
where /?* = (acr -1 , —a 2 f2r -2 ,0) and a = c are the 
space and time shifts of the acoustic metric. By set¬ 
ting 4- = tpi(r,t)e im ^e ikz , II = tti {r,t)e im ^e ikz , Si = 
{r,t)e lm ^e lkz , S 2 = im’L, and S 3 = ik^, where ( k,m ) 
are the axial and azimuthal wavenumbers, the hyperbolic 
system reads as follows: 

<9(7Ti + cd r (Ci — airi/r ) = (ac — ima 2 12) n\/r 2 + 

+c ( k 2 + m 2 /r 2 ) ipi — cCi/r 

dtipi ~ cd r (ai/q/r) = (ac — ima 2 fl) ipi/r 2 — ck\ (5) 
<9 t £i + cd r (tti — aCi/r) = 2ima 2 Shp\/r 3 — ima 2 Q^i/r 2 

The set of Eqs. © is augmented with the constraint 
\C\ = \d r ^\ — Ci | = 0, which is used to monitor the qual¬ 
ity of the numerical results. One-way inward propagation 
from the horizon is accounted for by an ingoing-radiation 
boundary condition, imposed through an excision tech¬ 
nique. Details of the numerical procedure will be given 
in a forthcoming publication El- 

Following the standard prescription for scattering pro¬ 
cesses in Kerr black holes El , the initial condition is 
chosen as a Gaussian pulse centered at r = ro and mod¬ 
ulated by a monochromatic wave, 

it>i (r, 0) = Aexp[-(r - r 0 + ct) 2 /b 2 - ia{r - r 0 + ct)/c ]| t=0 . 

( 6 ) 

The corresponding power spectrum is a Gaussian distri¬ 
bution P(ui) = P max exp [— (to — cr) 2 6 2 /4c 2 ], centered at 
frequency a with spectral width 1/6. The superradiant 
regime is0<<r/f2<m for m ^ 1. 

As already remarked, the main purpose of our calcu¬ 
lations is to assess the amount of energy that can be 
extracted from a BEC-like sonic hole as a function of its 
angular velocity 12. We use as reference a set of param¬ 
eters relevant to a BEC of 87 Rb atoms Ell , with vortex 
core radius a ~ 0.2 //in and angular speed 12 ~ 18 KHz. 

It is interesting to note that the inverse transversal time 
of the BEC vortex, c/a ~ 15 KHz, is very close to the cor¬ 
responding value for a cosmic black hole of radius a ~ 10 
km. Such a quantitative match stems from the very low 
speed of sound in BEC’s, of the order of a few mm/s. 

In the following we take c = h/(y/2M£) and a = £, 
where £ = (87r pa s ) 1 ^ 2 is the healing length. The inte¬ 
gration of Eqs. © is performed in the space-time do¬ 
main r £ [0, 150] and t £ [0, 150], in units of a = 1 and 
a/c = 1 in space and time. The angular frequency is 
analysed in the range 0.14 <12 a/c < 14, corresponding 
to a frequency range 1.8 < 12 < 18 KHz and a density 
range 5 x 10 13 < po < 5 x 10 14 cm -3 . The initial Gaus¬ 
sian pulse is centered at ro = 50a and a = 0.512, with 
amplitude A = 0.3c and variance 6 = 10a. We perform 
our study for m = 1 and k = 0.02/a, corresponding to a 
condensate axial extent H = 0.9 mm. Violations of the 














Figure 1: A density map of the real part of i/q in the r—t plane 
(for the sake of clarity, the case of a quasi-localised wavepacket 
at m = 0 is being shown). Due to the curved background, 
the wavepacket trajectory bends towards the sonic horizon at 
r = a. The light cone becomes parallel to the r = a axis, 
since no signal can escape from the horizon. 


Figure 2: Time evolution of the energy gain (top) and cor¬ 
responding fluxes (bottom) of the wavepacket for a = 0.7c/a 
and = 1.4c/a. The solid lines correspond to the superradi¬ 
ant case (m = 1) and the dashed lines to the non-superradiant 
one (m = 0). 


constraint C(t) =0 are monitored over the entire space 
domain and are found to be consistently below 10~ 6 for 
all sets of parameters under investigation. 

Fig^shows a density map of the real part of ijji in the 
range r/a £ [1, 20] and tc/a € [0, 10]. The initial Gaus¬ 
sian pulse moves towards the vortex horizon placed at 
r = a and its trajectory is bent by the potential outside 
the horizon. The bending of the trajectory, with light- 
cones heading towards the horizon, is consistent with sim¬ 
ilar findings in numerical relativity nm as shown in the 
diagram in FigH| 

In Fig0 we show a typical time evolu¬ 
tion of the energy of wavepacket E p (t) = 

(poM/2) / 0 27r dc/) dz f^ 43a v\r dr with v\ = 

(top curves), normalized to its initial value E p ( 0), as 
well as the (independently calculated) rate F(t) of 
change of the energy (bottom curves), normalized to 
its initial value F( 0), for a = 0.7 c/a and fl = 1.4c/a, 
within the superradiant regime (to = 1 , solid lines) and 
outside it (to = 0, dashed lines). F(t) includes the net 
flux across the surfaces at r = a and r = 143a as well as 
a term due to the bulk compressibility, 
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In the non-superradiant case, the energy of the scattered 
wavepacket goes asymptotically to zero, indicating that 
all the energy of the impinging wavepacket is lost to the 
vortex sink. In the superradiant case instead, the energy 
of the back-scattered wavepacket exceeds its initial value, 
indicating extraction of energy from the ergosphere at the 
expense of the rotational energy of the vortex. Consis¬ 
tently with this picture, the energy flux for the superra¬ 
diant (non-superradiant) case lies above (below) its back¬ 
ground value during the scattering event, approximately 



Figure 3: Time evolution of the energy gain for different val¬ 
ues of Q all in the superradiant regime. Here m = 1, b = 10a, 
ro = 50a, and a = 0.5D. 


in the range 35 < tc/a < 55. The energy gained via su¬ 
perradiance is by no means small, as it is seen to exceed 
in this case twenty percent of the initial value E p ( 0). 

It is now of great interest to examine the dependence 
of the superradiant energy gain on the angular speed of 
the vortex, so as to possibly identify an optimal value at 
which such energy gain can be maximised. In Fig0 we 
show the time evolution of the energy gain E p {t)/E p { 0) 
for a series of values of f 1 in the superradiant range 
Qa/c £ [0.6, 14], as can be experimentally achieved by 
varying the density of the condensate. A sharp increase 
of the energy gain is observed in the region fl > c/a. This 
is plausible, since > c/a marks a transition from the 
regime where the radius of the ergosphere remains within 
a factor two of the sonic horizon, to the regime where it 
grows linearly with II, thereby creating a sizeable ergo- 
spheric shell where energy can be extracted from. 

Since BEC vortices are quantized, one is naturally led 
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Figure 4: A_E P is plotted in units of poMc 2 a 2 H /2 on a log¬ 
arithmic scale as a function of Sla/c for three values of a/SI 
(solid lines) and compared with the background vortex energy 
Eh on the same scale and in the same units (dashed line). 


to ask whether the wavepacket may annihilate the vortex 
by extracting all of its energy in a much more demand¬ 
ing ’break-even’ condition, that is E p {oo) — E p (0) = E\, 
where Eb denotes the energy of the background vortex. 
In Fig. Q]we show the background energy E\> (dashed 
line) and the total energy gain A E p = (E p { oo) — E p (0)) 
for three values of a/SI in the superradiant range(solid 
lines), as functions of Cl. In the perturbative regime (for 
Q < 3c/a, say) the efficiency of energy extraction from 
the vortex grows much faster with f l than the quadratic 
increase of the background energy. This is especially true 
at large values of the ratio a/SI. Although substantial 
values of AE p /E^ are - by definition - beyond the scope of 
the perturbative Klein-Gordon description used through¬ 


out this work, it appears that nonlinear effects may pri¬ 
marily determine the way in which the energy extraction 
behaves as it becomes comparable to the background en¬ 
ergy. The possibility that substantial superradiance effi¬ 
ciencies, as they emerge from the Klein-Gordon analysis, 
may persist even in the non-perturbative quantum regime 
described by the GPE cannot be ruled out. Even though 
the condition AE p = E^ may remain out of reach for a 
single wavepacket, one may still conjecture that a train 
of wavepackets could reach the goal. It would be interest¬ 
ing to test this conjecture by numerical and experimental 
means. 

As a further development of the present model, it will 
also be interesting to consider a vortex with no drain m 
JUj and to apply our analysis to superradiant scattering 
from a giant vortex [2l|, yj, [2j|. Such vortices have been 
found to have up to 60 quanta of circulation and can 
therefore be well approximated within the classical limit. 

In summary, numerical simulations of the Klein- 
Gordon equation for sonic perturbations impinging on a 
BEC-like vortex suggest the possibility that, under typi¬ 
cal conditions of BEC experiments, a significant fraction 
of the vortex energy may be extracted via the mecha¬ 
nism of superradiance. It would be very interesting to 
test the realisability of such an exciting scenario both 
via the numerical solution of the Gross-Pitaevski equa¬ 
tion and by actual experiments on rotating Bose-Einstein 
condensates. 
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